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Abstract 

Localised defect modes generated by a finite line defect composed of several masses, em- 
bedded an infinite square cell lattice, are analysed using the linear superposition of Green's 
function for a single mass defect. Several representations of the lattice Green's function are 
presented and discussed. The problem is reduced to an eigenvalue system and the properties 
of the corresponding matrix are examined in detail to yield information regarding the number 
of symmetric and skew- symmetric modes. Asymptotic expansions in the far field, associated 
with long wavelength homogenisation are presented. Asymptotic expressions for Green's func- 
tion in the vicinity of the band edge are also discussed. Several examples are presented where 
eigenfrequencies linked to this system and the corresponding eigenmodes are computed for 
various defects and compared with the asymptotic expansions. The case of an infinite defect is 
also considered and an explicit dispersion relation is obtained. For the case when the number 
of masses within the line defect is large, it is shown that the range of the eigenfrequencies can 
be predicted using the dispersion diagram for the infinite chain. 

1 Introduction 

Despite being first studied in the late 17 th century by Newton pLj, wave propagation through 
discrete structures remains an active area of research today. A well known and interesting fea- 
ture of discrete media is the existence of pass and stop bands. The present paper examines the 
effect of a finite line of defects in an infinite square lattice. The behaviour of a lattice with a sin- 
gle point defect, or point source, can be described by the lattice Green's function. Such Green's 
functions have been studied in [2] for the two-dimensional square lattice. The resulting solu- 
tion was analysed for frequencies within the pass band and the corresponding asymptotics at 
infinity were also obtained. 

Movchan and Slepyan 1 3 1 examined several classes of continuous and discrete models with 
various forcing or defect configurations. When the forcing frequency, or natural frequency of 
the defect, is located in the stop band, localised modes were identified. For a particular choice 
of the mass variation, these defect modes can then be linked to stop-band Green's kernel which 
can be used in the construction of the defect modes as discussed. 

Lattice Green's functions are often studied in isolation and have proved a rich area of re- 
search (see, for example, 0], f5], Qjj, and references therein). For d -dimensional lattices, the 
Green's function is typically expressed as a d- dimensional Fourier integral. It is often possible 
to evaluate one or more of the integrals, as in f3] , but for d > 1 the Green's function cannot be 
expressed in terms of elementary functions. In the present article, several different represen- 
tations of the square lattice Green's function are presented, which prove useful for band edge 
expansions. 

Classical applications in the theory of defects in crystals and dislocations follow from the 
fundamental work of Maradudin 1 7 1 , where explicit closed form solutions were derived for a 
heterogeneous lattice system when two distant particles of different masses are interchanged. 



More recently, the envelope function based perturbation approach was developed in [8] and [9J 
for analysis of waveguides in photonic crystal structures. In the latter case, an array of cylinders 
(inclusions) represents a waveguide within a two-dimensional structure, and the frequencies 
of the guided modes are close to the band edge of the unperturbed doubly periodic system. 

Localisation of waves due to an infinite line defect embedded in an infinite square lattice, 
has been considered in the paper by Osharovich et al. [10] • For the case of an infinite line de- 
fect, dispersion relations can be computed in explicit form allowing spatially localised waveg- 
uide modes to be analysed. 

The book by Slepyan [ 1 1 J presents a detailed discussion of applications for dynamic lat- 
tice problems involving cracks modelled as semi-infinite faults, for both square and triangular 
elastic lattices. For a structured interface and a crack propagating with constant speed within 
a square lattice, localised modes were analysed in [12]. In particular, it was shown that the 
crack propagation can be supported by a sinusoidal wave localised along the crack, which the 
authors refer to as a knife wave. Using the lattice model, the dispersion relations for the crack 
within the square lattice can be derived. As shown in numerical experiments, these relations 
allow for the prediction of the average crack speed within the lattice when a fracture criterion 
forthe crack path bonds is introduced. More recently, Nieves etal. [13] studied the propagation 
of a semi-infinite dynamic crack in a nonuniform elastic lattice. The crack stability was anal- 
ysed and it was shown that information regarding unstable crack growth could be obtained 
from the study of the steady state regime. 

For the finite-frequency regime, a theory of asymptotic homogenisation for square cell lat- 
tices has been implement by Craster et al.[l4\. This theory makes use of information related 
to standing wave modes found in the lattice problem. Then a two-scale asymptotic procedure 
can be applied in order to obtain an effective partial differential equation for the correspond- 
ing macroscale that contains information about the microscale structure. 

The paper [15J showed that point forces acting at saddle-point frequencies within square 
and triangular lattices produce localised primitive wave forms within the lattice. Similar lo- 
calised primitive waveforms were demonstrated for the in-plane motion of elastic lattices [16] . 
The shape of these waveforms, created by a point force, were linked to the dispersive properties 
of Bloch waves in the lattice. 

The structure of this paper is as follows. In section[2] the problem of a finite line of defects 
(created by a perturbation of point masses) embedded in an infinite square lattice is consid- 
ered. Several representations for Green's matrix are presented, including integral forms and 
representation in terms of a generalised hypergeometric function. Localised defect modes for 
the finite line are analysed in section [2J] Therein, the necessary and sufficient condition for 
the existence of localised modes is formulated, and asymptotic expansions in the far field are 
also presented. Band edge expansions are constructed using an analytic continuation of the 
Green's function. Illustrative examples for a finite number of defects are given in section [3] 
where eigenfrequencies and eigenmodes are presented and compared with asymptotic results 
from the previous section. For the finite line defect it is observed that, in contrast to the ID 
and 3D cases, a localised defect mode may be initiated by removing any amount of mass from 
a line of nodes in the lattice. The analysis of a finite sized defect region is accompanied by the 
waveguide modes that may exist in a lattice containing an infinite chain of point masses, as 
in section|ij The governing equations for such a waveguide and the solvability of the problem 
are discussed in section|4j In section 4.3 the dispersion relation corresponding to the localised 



mode for the infinite chain is given. Finally, in section|5j a numerical simulation illustrates that 
the solution for the problem of the infinite chain can be used to predict the range of eigenfre- 
quencies of localised modes for a finite but sufficiently long array of masses representing a 
rectilinear defect in a square lattice. 
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2 A finite inclusion in an infinite square lattice 

Consider a square meshing of M 2 such that each node is labelled by the double index neZ 2 . 
Let there be N defects (with N e No) distributed along n 2 = as shown in figure[TJ The defects 
are characterised by a non-dimensional mass < r < 1, where the mass of the ambient nodes 
is taken as a natural unit. The stiffness and lengths of the lattice bonds are uniform and taken 
as further natural units. All physical quantities, such as the frequency and displacement, have 
been normalized according to these natural units and are therefore dimensionless. 
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Figure 1: A finite line of defects in an infinite 
square lattice. The length of the links, the stiffness 
of the bonds and the mass of the black nodes is 
taken as natural units. 



Let u„ denote the complex amplitude of the time-harmonic out-of-plane displacement of 
node n. Then, the equations of motion are 



jV-1 



= (l-r)w 2 5 0> n 2 ^ u n 5 Ptni , 

p=0 



(1) 



where <x> is the radian frequency, e\ = o>2,i] T > and 5tj is the Kronecker Delta. By means of 
the discrete Fourier Transform 



neZ 2 



(2) 



the governing equation (TJ may be written 



N-l 



(<w 2 -4 + 2cos£i + 2cos£ 2 )w FF (f)= ^ u pfi e~ ip ^ 

p=0 



(3) 



The positive root of the parenthesised term represents the dispersion equation for the ambi- 
ent lattice. It is observed that for to 2 > 8 there exist no real solutions to the dispersion equa- 
tion. Hence, the ambient lattice possesses a semi-infinite stop band: <x> 2 e (8,00). Inverting the 
transform yields the discrete field 



N-l 



w„M = (l-r)<w 2 ^ u Pi0 g{n,p;co), 

p=0 

where g{n, p; co) is the shifted Green's matrix defined as: 



(4) 



iff cos([ni-p]^i)cos(n 2 ^2) , r , e 

g{n,p;<o)=^ — — — d£id£ 2 - 

7r z J J co z -4 + 2cosd;i +2 cos 4 2 





(5) 
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For the purposes of numerical evaluation and asymptotic analysis in the stop band of the am- 
bient lattice {co 2 > 8), it is convenient to rewrite the Green's matrix as a single integral 



g{n,p;co) = 



1 

271 



(V^ 



1 



l)\ni-p\ 



it 



cos(n 2 £ 2 )d£ 2 , 



(6) 



1 



where a = co 2 /2 — 2 + cos E, 2 - Reversing the order of integration yields the same result, but with 
n\ — p and n 2 interchanged, and E,\ interchanged with E, 2 . An alternative representation can 
be found in the book by [17] , among others, as 



g{n,p;co) = 



{-1)" 



i-p+n 2 



I ni - p {x)I„ 2 {x)e ax dx, 



(7) 



where I m {x) is the modified Bessel function of the first kind, a = oo 2 /2 — 2 > 2. The integral 
is symmetric about ri\ — p = and n 2 = and therefore it may be assumed, without loss of 
generality, that ii\ > p and n 2 > 0. The integral J7) may then be represented in terms of a 
regularised generalised hypergeometric function (see [ 18|, section 3.15.6, equation 8) 



g{n,p;co) = 



l-iy 



: {{m + n 2 W 4 ¥ 3 



fli, a\, a 2 , a 2 
b\, b 2 , &1 + &2-I' 



(8) 



(2a) 1+m +"2 

where m = ri\ — p, a\ = (1 + m + n 2 )/2, a 2 = {2 + m + n 2 )/2, b\ = 1 + m, and b 2 = l + n 2 . The 
series |8) is convergent for a 2 > 4, that is, everywhere in the stop band of the ambient lattice. 
It is observed that along the ray m = n 2 , the Green's matrix may be written in terms of Gauss' 
hypergeometric function. In particular, equation |8) reduces to 



g{n,n,0;co) = 



((2/1)0 



112 



1/2 + n, 1/2+ n 
l + 2n 



(9) 



{2af+ 2n 

The function J9) is strictly positive in the region n > and a > 2. Hence, for a single defect, the 
lattice nodes along the diagonal rays do not oscillate relative to each other. 

Furthermore, for the case of m = n 2 = 0, the integral representation (5) reduces to the 2- 
fold Watson integral (cf. [4J and 1 6 1 among others) . Using a simple change of variables |5) may 
be written in terms of an elliptic integral, or alternatively, one may use J9) and observe that 

^1/2, 1/2 4 " 



g(0,0,0;w)=— 2 Fi 
2a 



= — K 
an 



a" 



(10) 



Together with equation (To) , the representation (7) is particularly useful since, by repeated 



integration by parts and use of the identity I n {x) = 21' Y {x] 
g(0, 0, 0; co) to a general g{n,p; co). 



7„_ 2 (x), one can iterate from 



2.1 Localised modes 



Of primary interest are localised modes, that is, modes of vibration at frequencies that are 
not supported in the ambient lattice and therefore decay rapidly away from the defect sites. 



Introducing the vector °U 
the eigenvalue problem 



[uo,o, u 2 fi,..., Un-i,oV and choosing n 2 = in equation |4) yields 



= {l-r)oj 2( S{co) a l/, 



(11) 



where the matrix entries [^(co)]ij = g{i - 1,0, j 
(and hence bisymmetric and centrosymmetric) 



l;co). Clearly, ^ is symmetric and Toeplitz 



^11 



#12 
^11 



^13 • 


• <Sl(N- 


1) 


'Sin 




■ Vw- 


2) 


^l(iV-l) 


^11 • 




3) 


^HiV-2) 




^11 




^12 



^11 J 



(12) 
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which greatly reduces the number of required computations. Indeed, for N defects the matrix 
'S has N independent elements. The solvability condition of the spectral problem i fTT) yields a 
transcendental equation in co, 



det \l N -{l-r)co 2( s] = 0, 



(13) 



where Ijv is the N x N identity matrix. Equation (13) is the necessary and sufficient condition 
for the existence of a localised mode. Symmetry implies that there exists an orthonormal set 
of N eigenvectors of 'S and hence, N eigenvalues (frequencies). The centrosymmetry of 'S 
allows the number of symmetric and skew- symmetric modes to be determined (see, for exam- 
ple, [ 1 9 1 ) . Introducing the N x N exchange matrix 



fo o o r\ 

10 

o ••' 

\J 0J 



(14) 



an eigenmode is said to be symmetric if = and skew- symmetric if = — J^ 6 !/ . For 
a system of N defects there exist [AT/21 symmetric modes and [N/2\ skew- symmetric modes, 
where [■] and [-J are the ceiling and floor operators respectively. Of course here, symmetry refers 
to the symmetry of the eigenmodes in the n i direction about the centre of the defect line. Due 
to the symmetry of the system, all modes are symmetric about the line n 2 = 0. 
Consider the total force on an inclusion containing N defects 



N-l 



F — (w p -i,o + Wp+i > o + 2u p ,ij . 



(15) 



p=0 



By definition, for a skew- symmetric mode u p ,o = — un-i- p ,o and further u Pt q = —ut]-\-p,q- 
Hence, for all skew- symmetric modes the inclusion is self-balanced (i.e. F = 0) and therefore, 
all skew- symmetric localised modes can be considered as multipole modes. 

For the illustrative examples presented later, the eigenvalue problem fTT] will be solved for 
the unit eigenvectors (|^| = 1). 



2.2 Asymptotics 

Here, asymptotics are considered for some particular cases. Asymptotic expansions for an 
isolated Green's matrix in various configurations have been considered in [3] and the approach 
detailed therein is used here. 



In the far field, along the line of defects. The case of n x — > oo, n 2 = and finite N is 
considered. Introducing the small parameter e = p/n\, the kernel of |6) may be expanded for 
small |e| <C 1. In particular, 



a 2 — 1 - ~ {^J a 2 — 1 - 1 — e log {^\J a 2 - 1 - a j 



il»il 



(16) 



it is observed that at large n\ and sufficiently small N, the dominant contribution to the inte- 
gral JB) comes from a small region in the vicinity of E, = %. Therefore, 



i 2 -l- 



j'" 1 P '~(Vc 2 -l-c) 



l»il 



\ni\ 



1 — s log yc 2 — l — c\ +s 



2Vc 2 -l. 



1 Iml 



(17) 
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where c = a> 2 /2 — 3. Thus, 



(V^-l-a) 1 " 1 ^ 1 

~ {^\J c 2 — 1 - cj 



lm-pl 



exp 



2\Jc z -\_ 



(18) 



In addition, 1 / \J a 2 — 1 ~ 1/ \J c 2 — \. Hence, for < £ <C 1 and making use of (6) 

,l«i-pl n 

>-£ 2 ) 2 



g(ni,0,p;w)' 



fVc 2 -l-cl 
27rVc 2 -l 



exp 



2\jc l -\ 



d? 2 



(19) 



Making the substitution x = {n — E, 2 ) \J\n\ — p\/2\J c 2 — X and performing the resulting integ 
tion yields 



g(ni,0,p;co)' 



(y/c*-l-c) 



\n\-p\ 



as ti] — > oo. 



^8n^c 2 -l V^l^i-Pi 
Thus from j4j, the physical field has the following approximate representation for «i — » oo 

, [Bl-Pl 



(20) 



^(V^l-c) 

w« 1 ,o(w)~(l-r)du / 2 J / / 

p=o \l%ns]c 2 -l Vl"i-Pl 



(21) 



where w p> o(a>) should be determined from fTT] . It is observed that when AT = 1 equation {2TJ is 
consistent with equation (4.17) of |3 up to a change in sign. 

In the far field, perpendicular to the line of defects. Here, the case considered is n x = 
p', n,2 — * oo with N and p' finite. The method used here follows the same general procedure 
as in the previous case. However in this case, the kernel is oscillatory and is therefore approxi- 
mated as a product of decaying and oscillatory functions. 

For sufficiently small \ p' — p\ and large n,2, the non-oscillatory part of the integrand in ([6) is 
approximated as before, leading to 



\"2\ 



g{p' ,n 2 ,p;co)' 



Uc 2 -\-c\ 

2n\j c 2 -\ 

n 

x exp 



(22) 



\n 2 \ 



2\[c l ~ 1 



cosQp'-plfOdfi. 



Making a similar change of variable, x = {n — £,\)\J \n 2 \/2\j c 2 — 1, and integrating, it is found 
that 



g{p',n2,p;co)~{-l) {p ' 



\n 2 \ 



c 2 -l 



x exp 



\n 2 \ 



(23) 



-ip'-pj 



\/c 2 -l 
2\n 2 \ 



Hence, for n 2 — » oo the physical field in (4) may be approximated by 



*V,« 2 M~(l-r)£o : 



(Vc 2 -l-c) 
\J%nsJ 



c 2 -l 



JV-1 



x J](-l/P'--exp 

p=0 



<p'-pY 



yV-l 
2\n 2 \ 



(24) 



u p ,o{co) 



\n 2 \ 
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It is observed that for N = 1 and p' = p, the above equation (24) is consistent with equation 
(4.17) of [3] up to a change in sign. Moreover, for the case ofp' = p (24) reduces to (21) . 



2.3 In the vicinity of the band edge 



The representations of Green's matrix (6) - (8) presented above are valid in the stop band. How- 
ever, given that the hypergeometric function in the representation (8) is zero balanced, that is, 
the sum of the bottom parameters minus the sum of the top parameters vanishes: 2{bi + b 2 ) — 
1 — 2(fli + a 2 ) = 0, the stop band Green's matrix can be extended to the boundary of the pass 
band by analytic continuation^ In particular, the analytical continuation of the function (8) 
has the form 



g{n,p;co) = 



( — 4) m + «2 

n{2a) l+m+n * 



{[l + m + n 2 ]/2)j 



x i 



if) 



jfc=o{([l + m + n 2 ]/2) ; } 



;${m,n 2 ,k) 



i/>(l+j-fc) 



,'l + m + n 2 .~\ , f 4 
+*P(l+j)-iP\ + H-iogl 1-^ 



(25) 



+(-iy'a)! X 



{([l + m + n 2 ]/2) fc } 



; $(m,n 2 ,k) 



where the reader is reminded that m = n\ — p, (•)/ is the Pochhammer symbol, tp{x) is the 
Digamma function, and 



~, {m) k {n) k 
${m,n,k) = 3F2 



(m + « 2 )/2, (m + « 2 )/2, -fc 
m, n 



;i 



(26) 



In this case, the continuation (25) holds for a 2 > 4, which in terms of frequency corresponds 
to oj 2 > 8. It is emphasised that in this section, the term "vicinity of the band edge" refers to a 
small interval 8 < a> 2 < 8 + e, where 0<s<l. 

Hence, choosing j = yields the leading order behaviour of (8) as a 2 — » 4+ {co 2 — > 8+), that 
is, as oj approaches the boundary of the pass band from the stop band: 



g{n,p;tL>)' 



n{2a) l+m+n * \ 



+ 



l+m + n 2 
2 



2y-xp 

1 

§{dl + m + n 2 ]/2) t } 2 



log 1 



(27) 



${m,n 2 ,k)\, 



where y is the Euler-Mascheroni constant. Alternative representations of the leading order 
continuations for general zero-balanced q +\V q were derived in (2T| . Since k > 0, the series rep- 
resentation of the hypergeometric function in (26) has a finite number of terms and therefore 
may be computed exactly. The convergence condition for the infinite sums in |25) and (27) is 
2+m+n 2 +j > 0, andis automatically satisfied since itwas assumed (withoutloss of generality) 
at the outset that m > and n 2 > 0. 

The asymptotic expression {27) is particularly interesting as it elucidates the nature of the 
singularity of the lattice Green's matrix at the band edge. In particular, the asymptotic repre- 
sentation (27) captures the logarithmic singularity as co 2 — > 8+. This logarithmically singular 
behaviour near the band edge is not obvious from the original representations presented ear- 
lier (cf. equations (6) -(H)). 

1 Indeed, for any integer balanced hypergeometric function q+iF q there exists an analytic continuation to the 
boundary of the unit disk (see |20|, among others, for details). 
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Figure 2: The solid red curve shows the 
asymptotic expression for the displace- 
ment field along the diagonal [n\ — 
n 2 with p — 0) in the vicinity of the 
band edge (cf. equation |28bfr ). The 
dashed blue curve shows the corre- 
sponding asymptotic expression for the 
field along the bond line (cf. equa- 
tion (28a)). The frequency chosen was 
co = 2.829. 



For some particular cases, equation (27) reduces to the following simplified forms. Along 
the rays m = (i.e. n\ = p) or replacing n 2 by n\ — p, along n 2 = 0: 



g(p,n 2 ,p;to)' 



(-4) 



l + «2 



n{2a) 



l+n 2 



2y + ip 



l + n 2 



+ log 1 



and along the diagonal rays m = n 2 (i.e. n\ = n 2 + p\. 



g(ni,m,p;(o)' 



16* 



n{2a) 



l+2m 



2r + 4> I - + m I +logl l 



= g (bond) (n 2 ;w), (28a) 



= g (dia g'(m;w), (28b) 



where the reader is reminded that m = n\ — p. The Digamma function grows logarithmically 
asm->oo and the term 2f + i/)(l/2 + m) is strictly positive for m > 0. Therefore, for sufficiently 
small m the bracketed term in equation (28) is negative in the neighbourhood of a = 2. Hence, 
in the vicinity of the band edge, the stop band Green's matrix exhibits fundamentally different 
behaviour along the bond lines compared with the diagonal rays. In particular, along the bond 
lines the masses will oscillate out of phase, whereas for the diagonal ray lines the masses will 
oscillate in phase, as illustrated in figure|2] In the far field, equations (28) further reduce to 



g(p,n 2 ,p;o))^ 
g{m,m,p;to)' 



(-4) 



i+n 2 



n (2a) l+n z 
16 m 



2y + log 



(!) 



+ log| l--r 

or 



n{2a) 



l+2m 



2^ + logm +log 1 



a' 



as n 2 — * oo, 



as m — > oo. 



(29a) 



(29b) 



Using equations (4) and |28) the anti-plane displacement for a lattice with N defects has the 
following asymptotic representation in the vicinity of the band edge 



N-l 



i ni ,o(<o) ~ (1 - rW " P ,og (bond) (ni - p; to), as cd 2 8 4 



p=0 
N-l 



( Bl ,oM ~ (1 - r)oj 2 " P ,og (diag) ("i - p; co), as ^ 8+, 

p=0 



(30a) 



(30b) 



along the rays n 2 = and n 2 = n\-p respectively. 



3 Illustrative examples 

Several particular cases are considered here corresponding to relatively short defects with N e 
[1,3]. The solid red curves in figure|3]show the / th solution, r w ,/(tu), of the solvability condi- 
tion (13) for a line of N defects. The shaded region indicates the stop band {a> 2 > 8) of the 
ambient lattice. For frequencies in this region, waves in the ambient lattice will decay expo- 
nentially away from the defect or source. It is interesting to note that in contrast to the ID and 
3D cases (see for example, |7|) the image of rjv,iv(w), indicated by the red curves in[3] is (0, 1). 
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In other words, a localised defect mode can be initiated by creating a defect in the lattice by 
removing any amount of mass from one or more nodes. In ID and 3D lattices, there is some 
upper bound on the ratio of the mass of the defect to the ambient lattice such that a localised 
mode can be initiated. As r — > 1, that is, the lattice approaches a homogeneous lattice, the fre- 
quency of the localised mode approaches the band edge {co 2 — > 8+). It is also observed that for 
N > 1, the solid red curves intersect the band edge at several distinct values of r. This suggests 
that for a given number of defects, there exists a maximum value of r below which all possible 
localised eigenmodes may be initiated. Above this value of r it is only possible to initiate a sub- 
set of the possible eigenmodes with the lower frequency eigenmodes being filtered out. In all 
cases, the highest frequency eigenmode persists for all possible values of r on (0, 1). For fixed 



co, the solvability condition 13 for a system of N defects is a polynomial, of at most degree N, 



in r. Therefore, there exist no more than N solutions for a given frequency co. 

The dashed blue curves correspond to the problem of an isolated chain of N particles of 
non-dimensional mass r*, connected by springs to two nearest neighbours and surrounded by 
rigid foundations. For such a problem, the out-of-plane displacement of mass «eZ satisfies 



se[vo,vir" ,v N -iV = 0, 



(31) 



where the matrix i£ has elements 



[Sfhj =(r*co 2 - 4)5 ij + Si-ij + 6tj-i. 



(32) 



The dashed blue curves in figure ^represent the solutions ; (o>) of the solvability condition: 
detif = 0. It is observed that as co — > oo, the blue curves approach the red curves from below. 
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(a) A single defect {N = 1) 



(b) A pair of defects (N = 2) 



(c) A triplet of defects {N = 3) 



Figure 3: The solid red curves show the I th solution, rjv,, (o>), of the solvability condition (T3) for a system of N defects 
embedded in the square lattice. The shaded region {co 2 > 8) indicates the stop band of the ambient lattice. The 
dashed blue curves show the corresponding i* solution, ,•(&>), of the solvability condition for an isolated system of 
N defects (cf. equation (32)). 



3.1 A single defect 

For the case of a single defect located at the origin, the quantity 'S in fTT) is a scalar: 



, , 1 f 4 
<g(o>) = —K[ -z 
an \ er 



(33) 



where K{x) is the complete elliptical integral of the first kind. The solvability condition maybe 
written as 



n,i = l + 7r 



CO 



1 

2 ~2 



K 



16 

(w 2 -4) 2 



(34) 
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which has the leading order asymptotic representation 



CO 



2 ' 



as 



CO — » 00. 



(35) 



It is observed that the solvability condition for equation (31) with N = 1 agrees precisely with 
the leading order high frequency asymptotic expansion. Hence, the coalescence of the red and 



blue curves in figure 3(a) 





Figure 4: (a) The localised defect mode for a single defect with r = 0.8 and co — 2.83. (b) The solid blue curve is 
the out-of-plane displacement along the line n 2 — and the red dashed curve is the asymptotic expansion for ri \ —> 
oo (cf. (2"T)). (c) The out-of-plane displacement along the line ri\ — n 2 (solid blue curve) with the corresponding 
asymptotic expansion (28b| for the band edge (dashed red curve), (d) As for (b), but the dashed red curve represents 
the band edge expansion along n 2 — (cf. equation (30a)). 



The localised defect mode is shown in figure 4(a) together with field along the line «2 = 
and the associated asymptotic field as n\ — > oo in figure 4(b) Figures 4(c) and 4(d) show the 
field (solid blue line) and the band edge asymptotics (dashed red) for a value of a = 2.006. The 
asymptotic expansions show good agreement with the computed field, even for the far field 
asymptotics in the neighbourhood of the defect. 



3.2 A pair of defects 

In the case of a pair of defects, ( §{to) is a 2 x 2 matrix with the diagonal elements given by (33). 
The off-diagonal elements have the form 



m W)]l2 = l --±K\^ 2 \. (36) 
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The solutions of the solvability condition are 



An(w 2 -A) 

r 2 ,i = 1 j jg ^, (37) 

w 2 n(w 2 -A)-2aj 2 {w 2 -8)K\ — - 7 -. 

\{w 1 — Ay 

Ant co 2 -A) 

ri,2 = 1 + j ^ r> (38) 

nco 2 {co 2 -A)-2(o 4 K\ 



Sco 2 -A) 2 

whence the leading order high frequency asymptotic expansions are 

3 , 5 

T% i ~ — ~ and Y2 2 ~ — o as a> — > oo, (39) 

which again, agree precisely with the solvability condition of isolated system (32) for N = 2. 
Hence, the coalescence of the red and blue curves in figure[3j 

Figure [5] shows the two defect modes together with the field along the lines n\ = 0, and 
n,2 = and the associated asymptotic field at infinity In addition, the fuchsia dash-dot line 



in figure 5(c) shows the band edge expansion in the vicinity of a = 2. In this case, figure 5(c) 



corresponds to value of a ss 2.025. Once again, the asymptotics are in good agreement with the 



computed field. Due to the symmetry, the field along the line n\ = 1 is identical to figure 5(e) 
for the symmetric case and identical up to a reflection in the line uo,n 2 = in figure 5(f) for the 
skew- symmetric case. 



The lower red curve in figure |3"(b) corresponds to rz,\ as defined in ( [37) . The maximum 



value of the lower red curve is given by 

r^= lim r 2 , 1 = 1 -. (40) 

oj^Vs 2 

Hence for a pair of defects, a symmetric localised mode cannot be initiated for r > 1/2. 

3.3 A triplet of defects 

For the case of three defects, the 3x3 matrix <^(co) has the [#]n and YS\ u elements as defined 
in equations ( [33} and J36). The remaining independent component is 

mw)] 13 = Mu)]n-^ + -E\X)> (41) 
2 7i ya J 

where E{x) is the complete Elliptic Integral of the second kind. The solutions of the solvability 
condition are of similar form to the previous two cases and are omitted for brevity. The high 
frequency asymptotics for r{to) are 

A-V2 4 4+v"2 
r3,i ~ r~> r 3,2 ?> and r 3>3 as <y->oo, (42) 

which again coincide with the solvability condition for ( [31] for the case of a particle triplet 

[N = 3). The maximum values of r 3 ,j(w) are r^™ = 1 — 37i/16, r^™ = 7/8 — (8 — An)" 1 , and 

(max) _ , 
'3,3 — 1 - 

For convenience, the three localised eigenmodes, along with plots of the associated asymp- 



totic expressions are shown in figures 10- 12 in appendix A. Plots of the displacement field 



along the lines n,2 = 0, n\ = 1 and n\ = are shown. The fuchsia dash-dot line in figure A.10(b) 
shows the band edge expansion in the vicinity of a = 2. In this case, figure A.10(b) corre- 
sponds to value of a ss 2.017. There are two symmetric modes (the lowest and highest fre- 
quency modes) and a single skew- symmetric mode, as expected from the properties of ^ dis- 
cussed in the previous subsection. However, for defects of mass r > r^^, it is not possible to 
initiate the lower frequency symmetric eigenmode and only a further symmetric mode and a 
skew- symmetric mode persist. For values of r > r^*\ it is only possible to initiate the highest 
frequency symmetric mode. 
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(a) Symmetric mode at co — 2.84 



(b) Skew- symmetric mode at a> — 3.35 





Figure 5: The localised defect mode for a pair of defects with r — 0.49. The blue curves are the out-of-plane displace- 
ment along the indicated line, and the red dashed curves are the associated asymptotic expansions in the far field (cf. 
equations (2T) and (24) as appropriate). The dot-dash fuchsia curve in figure 5(c) shows the band edge expansion (cf. 
equation (30a| ). 
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-6-5-4-3-2-1 1 2 3 4 5 n x 



2 
1 

n 2 = 
-1 
-2 
-3 
-4 
= 6 



Figure 6: A square cell lattice contain- 
ing an infinite chain of defects with non- 
dimensional mass r along n — 0, and an 
ambient lattice composed of particles with 
unit mass. As before, the stiffness and 
length of the links are taken as natural 
units. 



4 An infinite line defect embedded in a uniform lattice 

A numerical study of localised defect modes in a square lattice has been presented in [TO]. In 
this section novel results are presented for homogenisation approximations of long defects. 
The dispersion equations for an infinite line of defects in a uniform square lattice are derived 
and discussed in detail. 

The problem of an infinite line of defects embedded in a uniform square lattice, as shown 
in figure [6j is now considered. Given the symmetry about the line n 2 = 0, it is convenient to 
reduce the problem to a half-plane problem, which may be formulated as follows. 

4. 1 Equations of motion 

The equation of motion for a particle with n e Z x Z + is 

d 2 u n 



dt 2 

and for ri\ e Z, n 2 = is 

d 2 U„ 1>0 

r ^t 2 ~ 



Un+ei + U n - ei + U n + e2 + U n - e2 —4u n , (43a) 



: "«i+i,o+ u Hl -ifi + u nii i + u ni -\ —Au ni fl. (43b) 



Since the geometry in figure [6] is periodic in n\, and time-harmonic solutions are of primary 
interest, the solution u n is sought in the form 

u n = U n2 e i{nih -" t \ (44) 

where co > is the angular frequency, k g M is the non-dimensional Bloch parameter (nor- 
malised by the length of the lattice bonds) in the ti\ direction, and U n2 is the amplitude. As- 
suming the form (44), the equations of motion (43 a) and (43b) may be written 

f/„ 2+ i + [/„ 2 _i-2f2 1 (/c,iw)f/„ 2 = (n 2 >l), (45) 
t/i + t/-i-2fi r (/c,i6j)t/o = 0, (46) 



where 



n a {K,z)=l+2sm 2 {K/2)+ : ^ . (47) 



az 



2 



Solution outside the infinite line defect (| n 2 \ > 0). The solution for n 2 > 2 is then sought 
in the form 

U n2 =A n2 U lt |A|<1, (48) 

and for n 2 < —2, 

U n2 = -A n2 U 1 . (49) 
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The condition |A| < 1 imposes the localised displacement field about the chain of masses along 
n 2 = 0. The case |A| = 1 gives a field which propagates sinusoidally with constant amplitude 
away from n 2 = in the transverse direction. Together, equations (45) and (48) imply 

A 2 -2fii(>-,ko)A+l = 0. (50) 

The solution of (50) is 

{fii(fc, ico) - sign(fii(fc-, ico)) ^ fii(fc, ico) 2 - 1 for |f2i(jc,i<w)| > 1 , 
±1 forf2i(^,iw) = ±l, (51) 

Q,\(K,ia>)±i^ 1 — rti(K,ico) 2 for |f2i(?c,ia))| < 1 . 

It follows from (51) that if |ni(K,ko)| > 1, then |A| < 1. Further, if |fii(fc-,ko)| < 1, then |A| = 1. 



4.2 Skew- symmetric modes 

The solution for the problem when the skew-symmetry conditions are imposed along n 2 = 0, 
is zero. Indeed, for the case of a mode, skew- symmetric about n 2 = 0, the condition u nii n 2 = 
— u ni - n2 is imposed, and displacements are zero along the defect. Therefore, the solution of 
1 43a) which is zero along n 2 = and either decays or propagates with constant amplitude at 
infinity, is the trivial solution. 



4.3 Symmetric mode 

Dispersion relation of the chain at n 2 = for the symmetric mode. For the case when 
symmetry conditions are imposed about n 2 = 0, it will be shown that the dispersion relation 
for defect modes supported by the infinite line defect is given by 



co { ~\k) = 



r(2-r) 



l + 2sin 2 (K/2) 



+ 



V 7 ! + 4(1 - rf sin 2 (K/2)(l + sin 2 (K/2)) 



1 1 1/2 



(52) 



This dispersion relation is determined in two parts. First, the symmetry conditions are im- 
posed about the line n 2 = and a system is derived which links the displacements along the 
rows n 2 = and n 2 = 1. Then, the solvability of this system is considered for various cases of 
A, and (52) is deduced. 



The system for the displacements U and U\. For symmetric modes, the condition u null2 = 
u ni - n . z is imposed for n 2 > 0. In terms of the amplitude field the symmetry condition for n 2 = 1 
is U\ = U-\, whence equation (45) for n 2 = 1 and equation (46) give the system 

f/ + (A-2f2i(/c,iw))[/i=0, (53) 

Ui-n r {K,ico)U Q = 0, 

where (48) has already been used. It is convenient to introduce the following matrix notation 

S^ r {K,\co)u =0, (54) 

with 



f A-2f2i(/c,ko) 



1 



V 



l 



and u = [Ui,U ] T . 



-Q r {K,ico) J 



For non-trivial solutions Uj, j = 0, 1, of 54 , it is required that 



det(5^ r (/c,i£o)) = 0, 
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which leads to 

-(A 2 -2f2i(fo,ko)A)+l=0. (55) 



A 

Together with (50), equation (55) yields the equation 

A-fi r (K,ko) = 0. (56) 

Solutions of equation (56). In this part, it is shown that for < r < 1, there exist no solu- 



tions of the symmetric problem when |A| = 1, whereas for |A| < 1, the dispersion relation 1 52 
can be retrieved. For the solvability of (56), several scenarios are now discussed in detail: 



The case of\X\ = 1. Here it is proved that there exist no solutions of (56) for |A| = 1 and 
co>0. 



First consider the case when A = ±1. According to 1 51 , this corresponds to D.i[K,ico) = ±1 = 
A. Substitution of this into (56) , leads to 

Oi(fc, ico) — n r {K, ico) = . 

It follows that fii(/e, ico) = VL t {k, ico) if and only if r = 1. Thus for < r < 1, equation ([56) has 
no solutions for co > 0. Therefore, there are no solutions with constant amplitude for the case 
A = ±1 and < r < 1. The case of r = 1 corresponds to an intact lattice and one would expect 
constant amplitude (Bloch wave) solutions to be supported. 

The case of complex A remains. For A to be complex, the condition |f2i(fc,icu)| < 1 has to be 



satisfied, and then |A| = 1. Since r ^ 1, and tt r {K,ico) is real for co > and k e E, then (56) has 
no real solutions for co. 



The case of \X\ < 1. Together, equations (50) and (56) yield 

O r (K-, ico) 2 - 2f2 r (/c, iw)f2i(/c, ico) + 1 = , 
which by (47), is equivalent to 

- r(r - 2)co 4 + (1 + 2sin 2 {K/2))co 2 - 4sin 2 (^/2)(l + sin 2 (K72)) = . (57) 

Equation (57) is a biquadratic equation in terms of co. Now, according to (51), |A| < 1 occurs 
when |f2i(/c,io>)| > 1. Furthermore, for |A| < 1 equation (56), implies that |f2 r (^',iw)| < 1. More- 
over, the inequalities |f2 r (fc,i<x>)| < 1 and |fii(/c,ia>)| > 1, lead to 

r -l/2 w (l) W <0)< r -l/2 w (2) W (58) 

together with 

co {1 \k)>co or co>co [2 \k), (59) 

where 

a> (1) M = 2|sin(K-/2)|, and co [2 \k) = 2^1 + sin 2 ((c/2) . 



Either of the inequalities (59) implies that |A| < 1, and when one of these inequalities is taken 
with {58) , the solutions of (56) for co should satisfy these conditions. Equation (57) should be 
solved subject to conditions (58) and (59) in order to determine the dispersion equations. 
First note that for localised modes inequalities {58) and {59) yield either 

co i2 \K)<r- 1/2 co {2 \K), forallKeIR, 

which leads to r < 1, or 

r~ 1/2 ft) (1) (K)<a (1) (ir), forall>ceIR, 
which is never satisfied for any < r < 1. 
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Roots of the biquadratic equation {57) . The solutions of (57) are as follows: 

2 



w (±) W = j — [l + 2sin 2 (/c/2) 
I r(2 - r) L 



Vl + 4(r - l) 2 sin 2 (^/2)(l + sin 2 (;c/2)) 



1/2 



(61) 



for < r < 1. Here, it is shown that a> (+ ) is not a solution of (56), whereas a/ )(?c) is a solution 
of this equation. 



The function o>( ). Here it is proved that a>( ) satisfies ( 58 1 and (59) 2 , and is therefore a root 
of j56). 

Consider first, the function a>'~). Using Young's inequality, 

\f\ + 4(r - l) 2 sin 2 (?c/2)(l + sin 2 (/c/2)) < 1 + 2(r - l) 2 sin 2 (?c/2)(l + sin 2 (;c/2)) , 
and so from |6ll 



< 



r(2-r) 
4 

r(2-r) 



(1 + sin 2 (fc/2))(l + (r - l) 2 sin 2 (?c/2)) 



(l + sin 2 (^/2)Xl + (r-l) 2 ), 



for any /ceR. The last factor on the right hand side is positive and convex, and for < r < 1 is 
less than 2 — r. Therefore 



co { ~\k) < r- L '*a> w (K) , for all k ■ 



-1/2^(2) 



(62) 



Now, the function inside the radical of (<y' )) 2 (see 1 61 1), can be written as 

(r - l) 2 (2sin 2 (fc-/2)+ l) 2 + 1 - (r - l) 2 , 
and this implies that for < r < 1 and ceR, 



(oj { -\k)) 2 > -(1 + 2 sin 2 (K/2)) . 
r 



This then leads to 



co { ~\K)>r- L/z o) w {K), K 



1/2,^% 



(63) 



It remains to show co^ satisfies |59) 2 - Due to the equality 1 — (r — l) 2 = r(2 — r), it can be seen 
that this function is concave and less than 1 for r < 1. Then using this fact, co^ can also be 
estimated from below: 

a) { - ] {K) > co [2 \k) , for all k g R . (64) 



Then 1 62 , 1 63 and 1 64 show that &>( ' satisfies inequalities 1 58 and ( 59 1 2 and is therefore a 
solution of J56l . 



The function cq( + \ Now it is shown that <yM does not satisfy (58j and is consequently not 
a root of J56). Indeed, since the function in the radical of (wM) 2 is always positive 



■ sin 2 (fc/2) . 



r(2-r) 

for ireR. Then, since r < 1, it can be asserted that 

o) [+ \k)< r- llz co [l \K), foriceM. 
Therefore, o>W does not satisfy (58) and is not a solution of (56j . 
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Figure 7: The dispersion relation d/ - ', given 
in equation (52) , plotted as a function of 
o 0.5 1 1.5 2 2.5 3 3.5 4 the normalised Bloch parameter k/ti for r = 

k/ji [m _1 ] 0.05, 0.25, 0.5 and 0.75. 
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Analysis of the dispersion relation when |A| < 1 for various r. If a> satisfies |f2i(/c,io>)| > 
1 and |f2,-(fc,io>)| < 1 for all k e R and is a solution of (56) for < r < 1 then the dispersion 
relation for waves in the infinite wave guide is given by {52) . 
It is observed that when r = 0, co^ = 0. Further, 



-(l + 2sin 2 (/c/2))- 
r 



1 + 8 sin 2 (fc/2) + 8 sinV/2) 
l + 2sin 2 (/e/2) 



+ 0(r), 



where the second term on the right-hand side is bounded and negative and the first term dom- 
inates for r — » 0. Thus, co^~\k) — > oo as r — » 0. In figure[7] the dispersion relation 1 52 is plotted 
for several values of r. The in-phase standing wave solution, of the form (44) , is always given 
when k = and corresponds to the minima of the dispersion curves. The frequency of the 
in-phase standing wave is 



co = 



(65) 



r(2-r)' 

whereas for the out-of-phase solution, at k = n corresponding to the maxima of the dispersion 



curves, is 



CO = 



r(2-r) 



3+ x/l + 8(l- r) 2 



(66) 



5 From an infinite inclusion to a large finite defect: The 
case of large N 

In this section, the objective is to show that the range of eigenfrequencies for which localised 
eigenmodes exist for the model described in section |2j can be predicted using the model of 
an infinite chain of defects considered in section [i] A defect composed of N = 20 particles 
of non-dimensional mass r = 0.25 is embedded within an infinite square lattice. The eigen- 
frequencies of the finite defect were computed using the method described in section [2] and 
are shown as blue dashed-dotted, and red dashed, lines in figure [8] In this figure, the eigen- 
frequency co m i n = 3.0374 corresponds to an in-phase standing wave solution, whereas the fre- 
quency co max = 4.9344 represents the out-of-phase solution. The maximum and minimum 
eigenfrequencies are indicated by the red dashed lines in figure|8j 

Since N is large, it is useful to consider the model of an infinite chain embedded in a square 
lattice. Expressions (65) and (66) predict the values of the frequency co for which there exist 
such solutions. For the numerical values above, the in-phase solution occurs when k = and 
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k/ti [m 



Figure 8: The dispersion equation (52), for 
the infinite chain, plotted as a function of 
the normalised Bloch parameter, for r = 
0.25, represented by the solid black curve. 
Also shown are the blue dashed-dotted lines 
corresponding to the the eigenfrequencies 
computed for a finite defect containing N — 
20 masses. The dashed red lines correspond 
to w„ 



co = 3.0237 and the out-of-phase solution occurs when k = n and co = 4.9432. These values 
of the frequency are close to those encountered in the problem of the finite defect for N = 20. 
Moreover, all the eigenfreqencies computed for the finite defect lie within the passband for the 
infinite defect, as shown in figure [81 





(a) In-phase mode for co — 3.037 



(b) Out-of-phase mode for co — 4.934 



Figure 9: The solid blue lines are the eigenmodes for the maximum and minimum eigenfrequencies for a finite line 
containing 20 defects. The envelope function defined in (73) is plotted as the dashed red lines. 



Figure [9] shows the plot of the eigenmodes for the maximum and minimum eigenfrequen- 
cies computed for the line defect containing 20 masses. The maximum eigenfrequency co max 
corresponds to the out-of-phase mode, whereas the minimum eigenfrequency co m i n gives the 
in-phase mode. 



It is remarked that both, the field in figure 9(a) and the envelope of the field in figure 9(b) 



resemble the first eigenmode of an homogenised rectilinear inclusion. Using this motivation 
the difference operator 



®p (Op — (')p+ei (Op-ei (')p+e2 Hp - '2 ^(")p > 

is introduced. 

Making use of (4) it is found that 

I —^ + 1 Wni,o = (l-r)^ u pfi (@„ 1>0 + 6J 2 ) g{n l ,0,p;co), 



(67) 



(68) 
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where n has been restricted to {n : < n\ < N - 1, ri2 = 0}. Since the lattice Green's matrix is a 
difference kernel (i.e. depends on the difference \ n\ — p\), 



CO 



JV-l 



+ 1 u nit0 = (l-r)^ u p>0 (@ p ,o + oj 2 ) g{m,0,p;co), 



p=0 



whence, and recalling from (T) that (@„ + co 2 )g{n,p, co) = 5 ni>p 5 nZi o, it is found that 
(@„ + r&> 2 )w n = 0, for n e {« : < n x < N- 1, n 2 = 0}. 



(69) 



(70) 



It is observed that for a sufficiently large inclusion, the field above and below the inclusion 
behaves as u ni ,i = "m-i ^ Au nii o, with |A| < 1, in a similar manner to the infinite inclusion. 
Hence, using (70) together with the aforementioned approximation yields 



Un 1 +i,o + u„ l -i,o-2u n ,o+ [rco 2 - 2(1 - A)] W^.o^O, 



(71) 



for < «i < N — 1. The first three terms on the left hand side of (7T) correspond to the second 
order central difference operator. Hence, introducing the continuous variable r\ = n\ (where 
the reader is reminded that the length of the lattice links has been normalised to unity) equa- 
tion J7T) is written as 

■ d 2 



drj- 



+ rco-2{l-X) 



u{rf}fa0. 



(72) 



The form of equation (72) suggests that the homogenised system is analogous to a string on 
an elastic foundation, with the constant 2(1 — A) characterising the effective stiffness of the 
foundation. It is emphasised that |A| < 1 and as such, the stiffness of the elastic foundation is 
positive. 

Consider the problem of an infinite inclusion. According to equations {51) and (56) , the 
value of A corresponding to the lowest eigenmode is A = 1 + ra> 2 /2. For this value of A, the 
second order derivative vanishes according to equation (72) . Moreover, for the displacement 
at infinity to be finite, ^(17) must be constant for all 17 and the exact solution for the infinite 
inclusion is obtained. 

For the finite inclusion, it is observed that the displacements at the endpoints are small (cf. 



figure 9(a) \. Hence, for a simple estimate it suffices to impose u{0) = u[N — 1) = whence the 



solution to (72) is 



i{rj) = u sin ( v 7 ro) 2 - 2(1 - A)r? j , with A = 1 + - 



(73) 



where q is an odd number and uq an arbitrary scaling constant. The first eigenmode corre- 
sponds to A = —0.1396, which is close to the mean value of A obtained from the full numerical 
computation (A = —0.1426). The approximation (73) for A = —0.1396 is plotted in figure [9(a) 



as the dashed red line. The same approximation is used to produce the envelope function 



plotted as the dashed red lines in figure 9(b) One may observe that this, relatively simple, 
homogenised model predicts the envelope of the field very well. 



6 Concluding remarks 

A comparative analysis of two classes of problems has been presented: localised vibrations 
around a finite size defect created by a line of masses in a square lattice and an infinite waveg- 
uide represented by a chain of masses embedded in an ambient lattice. 

Although the physical configurations and the methods of analysis of these problems are 
different, one may observe remarkable properties of solutions, which can be used to make 
a strong connection. As illustrated in figure |8j the pass band for frequencies of waveguide 
modes, localised around an infinite chain of masses in a square lattice, contains all eigenmodes 
describing vibrations localised around a rectilinear defect built of a finite number of masses 
embedded into the lattice. 
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Special attention is given to the band edges: figure [8] shows that the frequencies of the 
eigenmodes for a finite rectilinear defect are distributed non-uniformly and they cluster around 
the edges of the pass band identified for the infinite waveguide problem. Furthermore, the 
limit, as one approaches the band edge frequency corresponds to a homogenisation approx- 
imation of the rectilinear defect as an inclusion embedded into a homogenised ambient sys- 
tem. The illustrative numerical simulation is produced for an array of 20 masses. However, the 
effect shown is generic, and, with an increased number of masses, the density of frequencies 
of localised modes near the band edges, identified for an infinite waveguide, increases. 

Symmetric and skew- symmetric modes have been constructed and analysed for a rectilin- 
ear "inclusion" built of a finite number of masses embedded into the lattice. It has also been 
shown that the total force produced by the vibrating discrete inclusion on the ambient lattice 
is zero for all skew- symmetric modes. Consequently, the displacement fields, associated with 
skew- symmetric modes, decay at infinity like dipoles, vanishing faster than the displacements 
corresponding to symmetric modes. This follows from the analytical representations for the 
solutions and illustrated in figures [5] and 10 where the skew- symmetric modes appear to be 



localised to a much higher degree than symmetric modes. In the aforementioned numerical 
simulations, the skew- symmetric and symmetric modes appear in pairs, and the frequency of 
the skew- symmetric mode is higher than the frequency of the corresponding symmetric mode. 
With reference to figure[3]it is also observed that, in contrast to ID and 3D cases, a defect mode 
can be initiated for any value of the contrast parameter r e (0, 1). In other words, removing any 
amount of mass for a point in a square lattice will yield a localised eigenmode. 

Finally, the reader's attention is drawn to the symmetric and skew- symmetric eigenmodes 
for a chain of 20 masses shown in figure |9j The corresponding frequencies are the maximum 
and minimum values among the array of frequencies associated with horizontal lines in fig- 
ure |8j The envelope curves for both diagrams in figure [9] represent the first eigenmode of a 
homogenised rectilinear inclusion. The simple homogenised model presented in section [5] 
provides the envelope curves for the finite inclusion. The form of the homogenised system 
suggests that, macroscopically, the inclusion behaves as a string on an elastic foundation. As 



expected, the skew- symmetric mode of figure 9(b) has the higher frequency than the symmet 



ric mode of figure 9(a) 
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A Appendix A. The localised field for a triplet of defects 

Following on from figurefTol the corresponding results for a triplet of defects is presented here. 



Figures 10 ■ 12 show the localised field for the case of N = 3 with contrast ratio of r = 0.4. In each 



case, the solid blue curves show the displacement field, whilst the dashed red curves show the 



associated asymptotics in the far field. The reader is referred to section 3.3 on pages 11 and 11 
for further discussion of the figures present here. 
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Figure 10: The first localised defect mode for a triplet of defects with r = 0.4. The blue curves are the out-of-plane 
displacement along the indicated line, and the red dashed curves are the associated asymptotic expansions in the 
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Figure 1 1: The second localised mode for a triplet of defects The blue curves are the out-of-plane displacement along 
the indicated line, and the red dashed curves are the associated asymptotic expansions in the far field (cf. equa- 
tions (2T) and (24) as appropriate). 
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Figure 12: The third localised defect mode for a triplet of defects. The blue curves are the out-of-plane displace- 
ment along the indicated line, and the red dashed curves are the associated asymptotic expansions in the far field (cf. 
equations (21) and (24) as appropriate). 
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